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Two-Dimensional Island Shape Determined by Detachment 
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Effect of an anisotropic detachment on a heteroepitaxial island shape is studied by means 

^1 , of a kinetic Monte Carlo simulation of a square lattice gas model. Only with molecular de- 

, position followed by surface diffusion, islands grow in a ramified dendritic shape, similar 

fS| . to DLA. Introduction of molecular detachment from edges makes islands compact. To un- 
derstand an anisotropic island shape observed in the experiment of pentacene growth on a 

• 1— I , 

O , hydrogen-terminated Si(lll) vicinal surface, asymmetry in detachment around the substrate 

, . 

I , step is assumed. Edge molecules detach more to the higher terrace than to the lower terrace. 

^ , The island edge from which molecules are easy to detach is smooth and the one hard to 

^ , detach is dendritic. If islands are close to each other, islands tend to align in a line, since 

\ detached molecules from the smooth edge of the right island are fed to the dendritic and fast 

^ , growing edge of the left island. 

' O ' KEYWORDS: heteroepitaxial growth, island shape anisotropy, detachment, lattice gas model, 

Q ' kinetic Monte Carlo simulation 

T— I ■ 

' 1. Introduction 

o : 

. For the production of high quality devices, it becomes more and more important to control 

■ atomic structures on clean crystalline surfaces. There are many studies, for example, on a shape 

' and density of two-dimensional (2D) islands which are nucleated and grown on a substrate 

\ surface by molecular beam epitaxy. They are mainly determined through the competition 

• . , between the diffusion of adsorbed molecules on a substrate surface and their incorporation 

■'^ ■ kinetics at the island edge. If the molecule incorporation at edges are very fast and diffusion 

\^ ' controls the growth, island edges undergo a morphological instability to form dendrites. In the 

■ " " ' extreme case of irreversible solidification, an island takes a ramified irregular shape with many 

branches,^ similar to a diffusion- limited aggregate (DLA).^' Its structure has no characteristic 
length and is fractal. On the other hand, if the incorporation kinetics is slow and limits the 
island growth, islands becomes compact. ^"^'^ 

On a vicinal surface steps exist and may bring about additional effects on surface struc- 
tures. For a homoepitaxial case, they lead to surface roughening via various instabilities as 
meandering and bunching. ^^"^^ For a heteroepitaxial case, substrate steps may affect a shape of 
epilayer islands. It is in fact observed recently in a heteroepitaxial growth of a pentacene (Pn) 
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submonolayer on a vicinal surface of hydrogen-terminated Si(lll) [H-Si(lll)] substrate. 
The experiment is performed in an expectation that vicinal steps give control on the epifilm 
orientation. Pn wets the substrate surface completely and form 2D flat islands in a submono- 
layer coverage. When a terrace width is large on a vicinal H-Si(lll) surface, Pn islands are 
isotropic and compact in shape. In contrast, when the slope of the vicinal substrate becomes 
steep and the step density increases such that an island covers many steps, the shape be- 
comes anisotropic; on the higher side of the vicinal surface an island is compact and rounded, 
whereas on the lower side it is dendritic. The island shape is found to be independent of the 
impinging orientation of the Pn molecular beam. As for the origin of this shape anisotropy, 
no definite conclusion has yet been stated. 

There is another experiment where anisotropy is observed in the shape of two-dimensional 
islands; a thin film of organic molecules on oxidized Ga-As substrate. The substrate has no 
steps, and the shape anisotropy is attributed to the tilting of long organic molecules in a 
crystalline film. Precise mechanism to determine shape anisotropy for this sytem still awaits 
for an analysis. 

Here we propose and study a possibility that may induce a shape anisotropy of islands; 
anisotropy in the rate of molecular detachment from an island edge. For example, on a vicinal 
surface, the rate of detachment from the edge of an island that climbs up a substrate step 
may differ from that from the edge of a descending island. Even on a flat substrate surface, if 
long molecules are adsorbed in a tilted conformation, detachment rate of molecules from the 
edge of an island may depend on a relative directions of tilting and migration. 

Whatever the mechanism may be, an anisotropy in the detachment rate should be reflected 
in the island shape. This is the aim we address ourselves in the present paper. We consider a 
simple and thus generic model where point molecules are adsorbing on a flat square substrate, 
and study it by kinetic Monte Carlo simulations.^'' The model is described in §2. In order to 
understand the effect of detachment on the island shape, we first study the case of an isotropic 
detachment in §3. In §4 we describe the possible origin of anisotropic detachment referring to 
the Pn experiment by way of example, and study the island shape in an anisotropic case. The 
last section §5 summarizes the results. 

2. Deposition and Diffusion on a Square Substrate 

Our model consists of point molecules depositing on a flat substrate surface of a square 
lattice. They diffuse on the substrate, and when they meet each other during the surface 
diffusion, they attach each other to form a dimer. When a diffusing molecule touches an 
already existing cluster, it attaches the cluster to increases its size, and the cluster grows to 
an island. When a molecule is deposited above an island, it diffuses on the island terrace 
till it reaches an edge, and steps it down to be incorporated into the island. Evaporation of 
adsorbed molecule from the surface is excluded, since the growth is performed under a low 
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temperature. With only surface diffusion and attachment at island edges, grown islands are 
ramified^ and look fractal, similar to DLA,^'^ at a low coverage. In order to realize compact 
islands, detachment or edge diffusion of molecules weakly bonded to the island edge are 
necessary unless evaporation is allowed. In the present model we allow detachment of those 
molecules connected to the island edge with only one nearest-neighboring bond, for simplicity. 
Still, they can be detached only to the substrate surface, but not allowed to hop up on the 
epitaxial island: There are three possible orientations for them to detach. 

We now review theories on an island density and shapes in terms of three parameters;^"^ 
a deposition flux a surface diffusion constant Dg, and a detachment rate Dg from a singly- 
connected edge site. As for the average island distance ^2 in a stationary state, we have to 
consider two-dimensional nucleation on a flat substrate. For a freshly deposited molecule, it 
takes about a time T2 = ^2/ to migrate a region of a linear dimension £2- During this time 
interval, other molecules are deposited to a density ni = Ft2, and the initial molecule has 
a probability ni^l to meet another molecule during its migration. (Logarithmic correction is 
necessary for more detailed analysis.) By assuming that the encounter causes the instantaneous 
nucleation, the encounter probability is the same with the nucleation probability per site. On 
the other hand, the latter is given by (a/^2)^ with a being the lattice constant, since there is 
one nucleus within a region of By equating the two results, one obtains ni^l = (0/^2)^ or 
a scaling relation 



This scaling holds when solidification takes place irreversibly and the critical size of a nucle- 
ation is ic = 1- In a more general case with a critical nucleation size i^, the island separation 
is given by^'^ 



where the lattice constant is set to unity. 

An irreversibly grown island with ic = 1 has a ramified profile similar to DLA, and consists 
of many edge molecules loosely connected only by a single bond. These loose molecules may 



J of breaking a bond. On a square lattice, the detached molecule can migrate three open 
directions on the substrate. If detachment happens along a straight island edge, it corresponds 
to an edge diffusion. If detachment takes place on the tip of a pointed branch, the detached 
molecule becomes a freely migrating molecule. But the fast terrace diffusion Dg may drive 
the molecule to another edge site not too far from the detached site. Therefore, detachment 
is effectively regarded similar to the edge diffusion. 

When the detachment is possible, the dendritic branches dissolve and the island can be 
compact in a polygonal shape. There is a theory on a stability limit of an edge length li at 




(1) 




(2) 



detach from the island with a rate = DgC "^/^^ 



at a temperature T with an energy penalty 
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which the polygonal shape becomes unstable.^' ^'^ On a straight edge of an island of a length 
ii, there should be a kink which incorporates molecules diffusing on the edge. With an edge 
diffusion constant D^, it takes about a time ri = if/Df, before a freshly deposited molecule on 
an edge is incorporated into the kink site. During this period, there are about FI'^ti molecules 
deposited within the territory of an island, where (.2 is an average separation between islands. 
With a fast surface diffusion D^, this number of molecules are deposited on the island edge. 
Thus, when it is of order unity, a one-dimensional nucleation takes place on the edge and the 
polygon side would become rough kinetically. Therefore, the critical size to keep an island in 
a polygonal shape is given as 




Fig. 1. An island shape at the coverage = 0.1 on a substrate with a size 500^. The ratio of the 
deposition flux F and the surface diffusion constant Dg is Ds/F = 10^^. The detachment rate 
D^/F is (a) 0, (b) 10^, (c) 10'', (d) 10^ (e) 10^ (f) 10^. A square at the left bottom is a maximum 
one that fits in the corresponding aggregate. 
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3. Isotropic Detachment 

In order to understand the effect of detachment on the island shape, we perform kinetic 
Monte Carlo simulations^'' of a system with deposition, diffusion and detachment processes. 
Before the study of an anisotropy effect, we start from the more fundamental case with an 
isotropic detachment in this section. To concentrate on the island shape variation, parameter 
values are so chosen as to allow only a single island in a large system with a linear dimension 
of -L = 500. In this large system, we provide initially a nucleus of a size 2 x 2 at the center of 
the system; This is the minimum size of an immobile island under a finite detachment rate, 
ic = 3. For a DLA-like island without detachment D^, = 0, the size of a critical island is = 1 
and the ratio of the surface diffusion Dg to the deposition flux F as high as Dg/F = 10^^ 
predicts a sufficiently large island separation as ^2 ~ 

{Ds/FY'^ ~ 300. A single island is, 
in fact, realized as shown in Fig. 1(a) at the coverage of = 0.1. The dimension df of the 
aggregate is estimated by assuming a scaling relation between the aggregate density c(r) and 
the distance r from the mass center as 

c(r) oc r'^f-'^. (4) 

If the aggregate is compact, the density is independent of the distance r and one obtains df = 
2. For a fractal DLA, its dimension is well-known to be df = 1.7, less than a space dimension 
2.^'^ Our simulation data is consistent with this value, though a deviation is observed around 
the center which is exposed for a long time under deposition; spaces between dendritic branches 
will be filled under a perpetual deposition, and the island becomes compact eventually. 

We now introduce detachment of molecules from the island edge to the substrate surface. 
The grown island at a coverage =0.1 becomes fat and anisotropic, as shown in Figs. 1(b) 
to (f). Every island consists of Na = 25,000 molecules. As D(,/F increases from 10^ to 10^ 
shown in Fig. 1(b) to (d), each branch becomes fat. On further increase of detachment rate 
to Df./ F = 10^ , the island becomes skeletal with a corner instability due to the Berg effect, 
as shown in Fig. 1(e). At a further large value of D^jF = 10^ (Fig- 1(f)) and more, the island 
remains square up to a coverage =0.1. 

On looking at island shapes in Fig. 1 one notices that there is a certain length scale which 
characterizes the width of dendritic branches, and this scale increases as the detachment 
rate increases. To quantify this characteristic length, one may think of various ways. The 
density c(r) is expected to show crossover from the central compact region to the outer fractal 
region. However, the central region is susceptible to small size effect and large fluctuation, 
and it is difficult to define a small crossover length definitely by this method. Another means 
is the box counting method used to characterize fractal objects. One covers the object by 
squares of a linear size b and counts the minimum number N(b) of boxes to completely cover 
up the object. It is proportional to for a compact object and to b^'^f for a fractal. If 
the object is compact up to a length scale ii, then the number N{b) is expected to show a 
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crossover behavior. In the present case, an island at Df, = shows a fractal behavior and that 
at Df./ F = 10^ a compact one, but for intermediate values of Df. the crossover turns out rather 
smooth and it is difficult to estimate the crossover length uniquely. Still another possibility is 
to investigate the edge length defined as the number of molecules with one to three nearest 
neighbor molecules. It is proportional to the island size Na for a fractal object and to ^JNa 
for a compact object. As the island size increases, one expects a crossover from compact to 
fractal. In fact, on the contrary, for small Df,/F < 10^, the perimeter length first increases 
in proportion to Na and then deviates from the linear relation. Because of the fast surface 
diffusion, islands first ramify, and then afterwards the profile is healed by a slow detachment 
process. 



10^ : 



10 : 




1 
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D,/F 

Fig. 2. Characteristic length £i of an island versus the detachment rate De- A slope represents the 

1 /2 

scaling behavior £i cx De . 

After all the incompetence, we estimate the characteristic length as a linear size of a 
maximum square box which can be completely covered by the island. Square boxes placed 
at the left bottom corner in Fig. 1 have the characteristic linear size ii thus determined. 
Due to the symmetry, there may be two kinks on an edge of a square island, and the length 
thus defined may corresponds to twice the characteristic length defined in eq. (3). In fact, we 
observe few kinks within the range of li on edges. But these kinks are transient, and we prefer 
the present way to estimate ii, since the main body of the island accumulates all the history 
of the growth and averages out fiuctuation. For the scaling law (3), in any case, a numerical 
factor is irrelevant, and we plot the length ii as a function of the detachment rate Df, in Fig. 2. 

At the coverage =0.1 in a system of a size = 500^, the island has a size Na = 25, 000, 
and the maximum linear size of a square island is \/Na ~ 158, which is indicated by a 
horizontal line in Fig. 2. For large Df,/F > 10^, the size £i reaches this maximum limit. In the 
intermediate range of 10"^ < De/F < 10^, the scaling behavior^ li oc dI^'^ of eq. (3) is found 
compatible with our data, as shown in Fig. 2. At small De, on the other hand, deviation from 
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the scaling law eq. (3) is apparent. This is related to our definition of £i. At small an island 
looks similar to DLA, and the outer dedritic branches in Figs. 1(a) and (b) are a very fine. 
However, the central portion of an island is exposed to a constant deposition flux all through 
the time until the coverage reaches to = 0.1. Therefore, the dendrite may thicken in the 
center. In contrast, for the DLA diffusing molecules start always far away from the aggregate 
itself, and central region is screened from the deposition flux. Thus, the central portion of the 
present island can cover-up large squares which are not covered in the outer region or in the 
true DLA. Since our estimation of £i requires only a single square box which is completely 
covered by an island, the present method overestimates ii, especially for small ii at small Dg. 
This causes the deviation from the scaling law at small D^. in Fig. 2. 

4. Anisotropic Detachment 

We now consider a shape anisotropy induced by an anisotropic detachment rate. The 
problem may be relevant to the shape anisotropy observed in the heteroepitaxial growth of 
Pn islands on a vicinal surface of H-Si(lll) substrate^^ or of an organic epifilm on a flat 
oxidized GaAs substrate. In the following a Pn system is chosen by way of an example to 
explain the possible origin of a detachment anisotropy, but the applicability of our generic 
model is by no means restricted to this specific system. 

Pn is a flat, elongated molecule, and on Il-Si(lll) surface it is known to grow a wetting 
film with the molecular long axis perpendicular to the surface, namely in a standing-up ori- 
entation.^^ It means that Pn interaction with the substrate surface is weaker than the Pn-Pn 
interaction. On a flat terrace, deposited Pn molecules make surface diffusion and they are 
incorporated into a 2D island when they encounter. Island are compact and their shape is 
isotropic. If the incorporation were instantaneous and irreversible, the island should be rami- 
fied and fractal, in contradiction to the experiment. A detachment from the island edge may 
accomplish an observed compact island shape on a flat terrace, as discussed in the previous 
section. We here study an asymmetry observed in island shape when Pn grows on a vicinal 
H-Si(lll) surface. As the step density in the vicinal surface increases and an island spans over 
many steps, island edges at the higher level are round whereas lower edges become dendritic. 

Steps may affect molecular diffusion on a substrate surface. If the step enhances an energy 
barrier against surface diffusion, the effect is symmetric around the step. The rate of step-up 
diffusion slows down by the same amount as that of step-down diffusion. Since no net flow 
is expected in this case, islands should be symmetric. If a step edge makes an extra bonding 
with a diffusing molecule, it will reside preferentially in front of the step. Then, nucleation 
rate will be enhanced in front of the step, and consequently islands will be nucleated along 
it, as often observed in various experiments. In the Pn experiment, however, this kind of 
island alignment along steps is not observed. Therefore, these kind of step effects on surface 
diffusion may be negligible for the Pn experiment. 
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Instead, we consider a step effect on molecular detachment from island edges. Imagine 
that a Pn crystalline film grows in a step-down direction and reaches an upper side of the 
descending step. Then, further crystal growth takes place by incorporating a Pn molecule 
from the lower terrace. In this case, the whole length of a newly attached molecule interacts 
with the island edge and the step ledge, and it is hard to be detached back to the lower 
terrace. The detachment rate is diminished by a factor Rl < 1. On the contrary, when the 
film grows in the step- up direction and reaches the lower side of the ascending step, the next 
Pn molecule to be incorporated is located on the higher terrace, and its upper part does not 
contribute to the Pn-Pn bonding. Since the substrate step has a height of about one fourth 
of the length of a Pn molecule, a loss in the bonding energy may be significant. Also, the free 
upper part of the long molecule may be susceptible to thermal fluctuation, and there may be 
a substantial entropy contribution which reduces the bonding free energy. Accordingly, those 
molecules incorporated into an island edge from the higher terrace are easily detached back 
with an enhancement factor Rh > 1. We now study the effect of this anisotropic detachment 
with Rh and Rl on the island morphology. 



(c) 



Fig. 3. Island shapes at different detachment rate from the left and the right edges. The surface 
diffusion is Dg/F = 10^^ and the detachment rate is D^/F ~ 10^. Detachment to the higher 
terrace (to the left) is enhanced by a factor Rh and to the lower terrace (to the right) is reduced 
by Rl, and theh values are {Rh,Rl)= (a) (2, 0.5), (b) (5, 0.2) and (c) (10, 0.1). The coverage is 
= 0.1 on a substrate with a size 500^. 

As we want to find out the essential factors of shape determination, the model here is 
chosen to be as simple as possible: Instead of providing a vicinal surface tilted down in +x 
direction, we simulate a model of heteroepitaxial growth on a fiat substrate surface, but with 
anisotropic detachment rate. The molecule that resides on the left edge of an island and 
connected to the island only by one bond has a higher detachment rate D^Rh with Rh > 1, 
whereas that on the right edge detaches from the island less frequently with the rate D^Rl 
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and Rl < 1. Those molecules residing on the upper and the lower edge of the island detach 
as before with the rate Df,. The diffusion constants are set to Ds/F = 10^^ De/F = 10^ to 
allow a single island in our system of a size = 500^. By initially providing an embryo of a 
size 2 by 2 at the center of the system, an island is nucleated and grows to a shape as shown 
in Fig. 3 at the coverage G = 0.1. The detachment enhancement factor to the higher terrace 
Rh is chosen to be inversely proportional to the suppression factor to the lower terrace Rl 
in order to decrease the number of free parameters, and they are varied as (a) {Rh, Rl)={'^, 
1/2), (b) (5, 1/5), and (c) (10, 1/10). In the isotropic case with Rh = Rl = 1, the island has 
a skeletal shape shown in Fig. 1(e). Even though enhancement and the suppression factors, 
Rh and Rl-, are applicable everywhere in the system and there is no local inhomogeneity, 
the left and the right halves of the island show a drastic contrast: The left-hand side (Ihs) is 
polygonal and the right-hand side (rhs) is dendritic. 

To quantify the shape difference in the Ihs and the rhs of an island, one may think to 
utilize fractal dimensions in both sides. In Fig. 3(c), for instance, the Ihs of an island is very 
narrow, and a box counting in a limited width yields only a rough evaluation as df ~ 1.9 ~ 2. 
On the rhs, dendrites are clearly observable but they are still fat and not long enough. A 
fractal dimension depends on the narrow fitting range as dj ~ 1.7 ~ 1.8. For a more accurate 
determination of dj larger systems are necessary, but it is out of scope of the present study. 

As for the upper and lower edges in Fig. 3, the Ihs and the rhs of an island also show 
a contrast. In fact, the detachment rates from them are the same Df, as in Fig. 1 (e). The 
Ihs of islands in Fig. 3 are similar to that in Fig. 1(e), whereas the upper and lower edges in 
the rhs of an island in Fig. 3(c) are quite flat compared to the rhs of an island in Fig. 1(e). 
It may be because branching of the upper or lower edges to the left ( or to —x direction) is 
suppressed. Branching takes place only to the right, and these dendrites grow to the right ( 
or +x direction ) fast. This suppression of left branching decreases fluctuation of upper and 
lower edges of an island. 

Furthermore, on looking through Fig. 3(a) to (c), one may get an impression that an island 
seems to be shifted to the right as Rl decreases. This apparent shift is due to the difference in 
the growth rate in ibx-directions. The dendritic right edge grows faster than the facetted left 
edge. Difference in the two edge speeds is quantified in a previous experiment^^ by a ratio of 
growth distances of side edges from the nucleation center. Since the initial nucleus is provided 
at the center of the system Xc = Uc = L/2, the degree of the shape anisotropy is measured 
by the maximum and the minimum distance in the x direction of the island edges from the 
initial nucleation center as 

q = Z • (5) 

o Xq Xmin 

The simulation result is summarized in Fig. 4. As the reduction factor Rl ol detachment to 
the lower terrace decreases, the shape anisotropy increases. 
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Fig. 4. Island shape anisotropy L/S versus detachment anisotropy Rl = 1/Rh obtained from Fig. 
3. Line is a guide for eyes. The surface diffusion is Dg/F = 10^^ and the detachment rate is 
De/F = 10^. 

So far, we discussed the shape of a single island. In a large area of the vicinal surface 
many islands are nucleated and they affect their growth and shape each other through the 
diffusion field. Their average separation, for instance, is determined by the ratio of the surface 
diffusion constant Dg and the deposition flux F, if detachment anisotropy is discarded, as 
described in §2. With anisotropic detachment, we obtained spatial distribution of islands as 
shown in Fig. 5. The surface diffusion constant in Fig. 5(a) and (c) is set to Ds/F = IQi^^ 
and in Fig. 5(b) to Dg/F = 10^. Detachment rate in Fig. 5(a) is set D^jF = 10^ with the 
enhancement factor Ku = 10 and reduction factor = 0.1. Islands are well separated and 
they seem to be distributed homogeneously. In Fig. 5(b), the surface diffusion is decreased 
to Ds/F = 10^, and the island density increases. Then one may detect some correlation in 
the alignment of islands in x-direction. In Fig. 5(c) the surface diffusion is kept the same 
with the case in Fig. 5(a) but, instead, the detachment rate at the upper and the lower edges 
of an island is increased by a factor of ten: Molecules are easily detached from the three 
edges except the right one, and accordingly the right dendritic edge grows faster. Enhanced 
molecule detachment in y direction enlarges the longitudinal island separation and weakens 
the correlation in this direction. In contrast, one has an impression that islands align laterally 
in X direction. 

In order to quantify structural correlation, structure function or the Fourier transform of 
the density correlation is usually employed. In the present small system, however, the structure 
function is noisy and it is difficult to identify the anisotropy of the island alignment. Instead, 
we use a simple measure to indicate an anisotropy in x and y directions; number densities 
of empty rows and columns which are free of adsorbed molecules. For example, in Fig. 5(a) 
4.8% of the vertical columns are free of adsorbed molecules whereas 9.2% of the lateral rows 
are free. Though each island has a clear anisotropic shape, islands arrangement in space is 
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rather isotropic. In Fig. 5(b), the empty row increases to the density of 16.6%, whereas there 
are no empty columns; there are at lease one adsorbed molecule on every column. In Fig. 5(c) 
alignment in x direction tends stronger; 51% of empty rows but no empty columns. Thus, the 
measure of empty rows and columns seems very clear and helpful to indicate the degree of 
anisotropic arrangement of islands. 





Fig. 5. Islands distribution at small surface diffusion, (a) and (c) Dg/F = 10^°, and (b) Ds/F = 10^. 
Other parameters are (a) and (b) De/F = IQ-^, Rh = 10, i?L = 0.1, and (c) DJF = 10^ Rh = 
1, Rl = 0.01. The system size is 500^ and the coverage is O = 0.1. 



5. Summary and Discussions 

An effect of an anisotropic detachment from island edges on the determination of a two- 
dimensional island shape is studied by means of a kinetic Monte Carlo simulation of a square 
lattice gas model of heteroepitaxial growth. The study may be relevant to the formation 
of anisotropic islands observed in heteroepitaxial growth experiments of Pn molecules on a 
vicinal surface^^ and/or of tilted organic molecules on a flat substrate. Only with a molecular 
deposition with flux F followed by a surface diffusion Ds, islands take a ramified dendritic 
shape, similar to DLA. Addition of detachment from the edge with a rate De lets dendritic 
branches fat, and with a large compact square islands are formed. 

By assuming that a detachment rate depends on an edge orientation, island shape becomes 
anisotropic. For example, the left edge with an enhanced detachment with a rate RhD^. with 
Rh > 1 is found smooth and compact, whereas the right edge with a small detachment rate 
RiDe with Rl < 1 is ramified and dendritic. If islands are close to each other, one notices 
a tendency that islands align linearly, since detached molecules from the smooth edge of the 
right island are fed to the dendritic and fast growing edge of the left island. 

We comment on the anisotropy factors, Rh and Rl, in Pn experiments.^^ Even without 
steps, Pn islands have compact shape, and thus the detachment process takes place often 
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enough. This means that a Pn-Pn bonding is not too strong to prevent edge molecules from 
detachment, and some part of the bonding is broken thermally. Since the Pn molecule is long, 
it is conceivable that the upper part is most susceptible to thermal fluctuation. When substrate 
steps are present, a Pn molecule incorporated from the higher terrace misses bonding at its 
upper part, but if thermal fluctuation already breaks the bonding there, the effect will be 
small. Then, Rh can be close to unity. For a Pn molecule attaching from the lower terrace, 
on the other hand, its upper part may make a strong bond with the pre-existing island edge 
which lies higher. Then, the detachment may be strongly supressed, yielding a very small 
Rl. Estimation of precise values of these factors requires microscopic information of chemical 
bondings, molecular rigidities, and so on, and out of scope of the present study. 

Above explanation is refering to the Pn epilayer on a vicinal surface just by way of example, 
but it is not meant to be specific to that system. We expect our model is universal and is 
applicable to other cases, such as heteroepitaxial growth of tilted molecules on a flat substrate 
surface. 
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